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1.  Formulation  of  the  problem  and  fundamental 

designations  *  “  , 

We  are  going  to  investigate  surface  chemical  reac-j 
; tions  in  solids ,  which  do  not  take  place  simultaneously 
everywhere  in  the  system,  but  begin  in  separate  points  of- 
:  :the  solid,  and  from  there  spread  throughout*  The  creatioh 
i of  “seeds”  of  "initiation  centers"  of  s  similar  reaction  ; 
:may,  as  a  function  of  conditions,  take  place  in  the  inter** 
!ior,  or  on  the  surface  of  the  solid,  on  the  edges,  ridgesj, 
■  or  angles  of  the  crystals* 

The  similar  surface  chemical  reactions  comprise: 
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•  1)  the  transformation  of  one  solid  phase  into  another  —  ; 

I  recryst’allitation,  and  2)  the  transformation  of  one  solid; 
i body  into  another  with  the  simultaneous  evolution  of  a 
I  gaseous2'product  j  as  is  evidenced,  for  example,  in  the  case 
i of  the  decomposition  of  carbonates. 

The  theoretical  analysis  of  the  kinetics  of  simil— j 
i ar' reactions  which  are  started  on  the  crystal  surface  has' 

I  encountered  considerable  difficulties  arising  from  the 
I  fact,  that  the  transformation  a ones  which  spread  from  the 
i various  •’initiation  c enters do  not  grow  freely,  but  nm 
i across  one  another  according  to  the  conditions  of  the  in- 
‘dividual  case.  In  the  studies  of  Bradley,  Colvin  and  . 
iHume  (1),  Roginsky  (2),  Ismailow  (3),  and  Iferof eyev  (4), 
ItMerefore,  only,  separate  and  individual  cases  were  invest- 
! igated  —  the  beginning  or  the  end  of  the  process. 

I  The  example  where  the  seeds  of  the  new  phase  are 

formed  in  the  icitex'ior  of  the  original  phase.-g,  accordir.g  ' 

I  to  the  requirements  of  the  case,  but  thereafter  grow  untif 
I mutual  saturation  was  recently  analyzed  by  Kolmogorov  (5)^ 
land  also  by  Johnson  and  Hehl  (6),  by  some  different  methods. 
'The  couclx'^sicns  of  these  authors  represent  some  modifica-' 
:tion  of  a  new  method,  which  as  far  as  we  know  was  first 
i  suggested  by  Clausius  during  his  .determination  of  the  laws 
iof  the  partition  of  the  free  path.. 

J  We  have  accepted  Clausius*  method  in  the  case  of 

i chemical  reactions  or  phase  transformations  initiated  on  i 
I  the  surface  of  the  solid.  In  this  work  we  shall  confine  ■ 

I  ourselves  to  the  case  of .  the  j’l&Jis  surface,  i.e.,  to  part- 
jicles  of  very  large  dimensions.  The  reactio.n  kinetics  j  - 
■for  particles  df  ordinary  dimensions,  whose  surface  has  ai 
; curvature  different  from  sero,  may  be  derived  in  a  similar 
iway.  ■  ^  . 

;  The  results  of  the  l.atter  investigation  will  be  pub¬ 

lished  by  us  separately. 

Jor  the  analysis  of  the  kinetics  of  the  process  for 
•;  infinitely  large  pa.rticles,  we  wrill  assume  that  from  every 
center  on  the  surface  there  spreads  into  the  interior  of  ' 
the  solid  a  spherical  front  of  the  process  with  a  constant 
linear  velocity  A  cm/sec,  Th©  initiation  centers  of  trans- 
tformation  may:  1)  already  at  the  very  beginning . of  the 
i process  be  found  spread  over  the  surface  of  the  body  at 
i ah  average  density  n  par  unit  area,  or  2y  these  centers  ! 
iare  created  as  a  function  of  time  with  an  initial  probabil¬ 
ity  a  per  unit  area  per  unit  time  interval.  Naturally, 

;a  more  complicated  case  is  also  possible,  such  that  the 
transformation  centers  appear  already  on  the  surface  part¬ 
ially  formed  as  a  function  of  time. 

At  the  beginning  of  the  process,  the  transformation 
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zones  are  uniformly  growingr,  hemisrneres.  As  time  fsoes 
on,  these  zones  be^^in  to  interpenetrate  one  another,  thus 
strongly  complicatirig  the  kinetics  of  the  process  and  its 
calculation.  Bventuaily,  after  a  sufficiently  long  period 
of  time,  all  the  separate  zones  merge  and  form  one  front 
of  propagation  of  transfcrniation,  which  travels  in  a 
direction  rerpendicular  to  the  surface  of  the  solid  with 
the  same  linear  velocity  . 

The  mathematical  a^nalysis  of  the  kinetics  of  the 
initial  period  —  the  period  of  free  zone  growth  —  does 
not  present  particular  difficulties.  jDeperiding  on  the 
nature  of  these  initiation  centers  (a  =  /  0,  or 

n^  =  0»  a  0)  we  obtain  the  yield  of  the  reaction  product 

proportional  to  the  third  or  the  fourth  power  of  time  t. 
The  kinetics  of  the  final  stage  is  simpler  still,  where 
after  the  formation  of  a  practically  plane  front,  the 
quantity  of  the  reaction  product  is  directly  proportional 
to  t.  The  analysis  of  the  intermediate  period  —  the 
interpenetration  of  zoiies  and  their  merging  into  one 
front  —  represer;ts  a  fairly  complicated  problem  in  the 
theory  of  probability,  which  as  yet  has  not  been  solved. 

2 .  Kine-iic  eoTiations  of  the  process. 

Considering  the  conversion  of  particles  of  suffi¬ 
ciently  large  dimensions,  we  say  neglect  the  surface  curv¬ 
ature  and  consider  them  as  planes  to  all  intents  and  pur¬ 
poses. 

Th\is,  in  the  limiting  case,  one  may  consider  the 
half-space  bounded  by  the  infirdte  plane,  and  filled  with 
an  isotropic,  homogerieous  substance,  capable  of  undergoing 
a  transformation. 

At  the  interface  v;e  have,  as  before,'  transformation 
centers  spread  on  the  surface  with  an  average  surface  den¬ 
sity  n^,  according  to  the  conditions  of  the  case,  and 

there  may  also  be  formed  new  cei:ters,  whose  probability 
of  formation,  taken  per  unit  area  and  unit  time  interval, 
we  adopt  as  a 

Let  us  consider  an  arbitrary  point  K,  v/hich  is  at 
a  distance  x  from  the  interface.  We  shall  define  the  pro¬ 
bability  of  the  point  K  finding  itself  in  a  transformation 
zone  at  time  t,  which  zone  is  propagated  from  some  trans- 
formation  centers  on  tixe  surface,  assuming  that  the  trans¬ 
formation  zone  spreads  from  the  center  spherically  with  a 
constant  linear  velocity  ^  .  naturally,  an  inequality 
t  ^  7i/Z  is  valid,  as  in  the  case  of  t  <dxAt  tlie  point  M, 
in  general,  will  not  be  in  the  reaction  zone. 


■  Thus  5  ws  will  tsks  li  Siid  .xn'/sst-it'a'ts  "the  t;iiae 

iKterval  (ti  t  dt  ■ .  Witbin  this_tiae  interx^al ,  prdy 
those  transformation  sorts  cia,';  ax-rdve  at  the  point  M  which, 
first?i.y,  propagate  from  tranafcr-uiation  centers  distributed 
with  a*^  surface  density  or?  the  circuififererice  of  a  circle 

of  radius  r  »  ~  'X~~  with  ring  of  wicith.  dr  (figure 

1);  and  secondly,  from  centers  foriaing  on  rings,  whose 
radii  satisfy  the  condition  ^  ~  yf  (figure  2), 

in  the  time  interval  dt,  which  is  counted*  ..nowev^g,  not 
from  time  t,  out  from  the  moment  ^  )* 

The  prohability  of  point  ¥■  fslling  iirco  the^  trans¬ 
formation  zone  during  the  time  interval  (t,  t  +  '.it/  Vir5.ll 
be  equal  to  the  sura  of  the  probability  of  falling  into  _ 
the  transformation  zone  propagating  from  centers  distri¬ 
buted  with  a  constant  density  n. 


p’  «  2/m  rdr  =  2/if'n_/^'^td’i 

‘  O  Q 


and  the  probability  of  falling  into  the  trsnsf cipration 
zone  pro'oagoting  from  centers  dui;lng  the  uliiie  intervsi 
dt,  but  a  little  earlier  is  tr 

•Y  t“  -  x-^ 

p"  =  dt*2/ra  *  \  r-’dr’  =:  rra(i"’'t‘"  -  x"^)dt. 


The  probability  oJ’  the  point  F  not  undergoing  t  ran '.'forma¬ 
tion  in  the  time  interval  (t,  t  -!■  at;  will  be 

Q  _  i  _  /  C'!  P^')  ^  1  -  ?ra(/l'^t"''  -  X 
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it'V  is  the  pi’obability  of 


‘t  +  -  x^lVdtj  . 


Expanding  the  left-harKi  si'de  into  a  pcwsr  series  in  dt , 
we  obtain  an  equation  determining  the  unknown  function 

f(x,t):  ^ 

».  ~/l/Tn  /ft  -  x^)  7dt. 

i(x,t ;  ~  o  - 


Fig .  1 


Fig.  2 


Solving  the  last  equation,  v;e  will  obtain  an  expression 
for  f(.,t):  -;ra(|A5  - 

x^t  +  |x^//)„7  » 

in  which  the  constant  oi‘  integration  is  found  from  the 
obvious  expression  f(x,  x/Z)  «  1.  We  have  obtained  the 
probability  of  the  point  M  not  undergoing  a  transforma¬ 
tion  at  cine  t.  Ihe  prob/oility  of  point  undergoing  a 
transformation  at  time  t  will  be 

p(x,t)  =  1  -  f(x,t)  =  1  -  exp/~/7hQ/?(t^  ~  - 

~  5tx-//^  +  .  (1) 

5  •  , , 

Let  us  distinguish  an  elementary  layer  of  height  dx  and 
with  base  area  of  1  in  which  point  n  is  contained. 

The  volume  of  the  substance  undergoing  transformation  in 
this  layer  at  time  t  is  equal  to  P(x,t)*l  cm^'dx,  and 
the  total  volume  experiencing  transformation  at  time  t, 
calculated  for  the  unit  surface  area  is: 

It  f  1 

v(t)  ^  J  ?(x,t)dx  =  7tJl  -  C  exp^rn^/?.^t^(l  -  2^)  - 


0 


0 


(1  -  +  2z7)  7<iz 


A 

J 


j 


(2) 


where  z-  =  x/Xt, 

In  the  px'eceding  discussion  the  quantity 
P(x,t)*l  cm^'dx  is  identified  with  the  portion  of  the 
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f  ■  ^ 


ime.  t  midepgoes  transformation  in  a 


fsnbstance  wMcJ^  at  tlL 
Tvolnma-  of  1  dsr*  / 

I  Tfee  transferms;^idn- rate  is 

s^)  ^ 

I  .(1  -  +  2z^)j^dz  I  ^  -  z^)  -V 

2  -.,  ■  ,  b 


•(1  -  5z  +  2z-’27exp^ffn^yi^t^(l  z^) 


•  (1.,.^-  +.2z^j7  dz*  ' 

I  avaiiiation  of  these  integrals  is  not  possible  ■ 
I  by.  conventional  means ^  and  it  may  be  done  approximately,  : 
only.  It  is  estpedient  to  analyze  first  the  initial  and  I 
.the  final  parts  of  the  process.  ; 


(D 


It 

.  » 


t 


Free  groxfth  of  separate  zones.  •. 

Til  the  initial  part  of  the  process,  i.e.,  for  staall 


t^l 


/ 


1 

'  I^o"? 


^  'r 

-  “n  •  ’  \ 


5/ 


ITaA"^ 


“  %g-  (ii:) 


lone  may  expand  the  function  under  the  integral  sign  in 
expression  (2),  liEiiting  it  to  the  first  three  terms. 

; This  approximation  yields: 

v(t)  ~  I 


cc 


'^'3  jjL./  ^ 


(1) 


give  the  sum! 

of  the  voluiaea  of  ail  srowth  zones  (hemisphere),  formed  ’ 
between  the  time  t  «  0  and  t  =  t,  tafeen  per  unit  area.  i 


The  quantities  and  ^jTa/ 


(t/tj^)'^  Slid  1  -  afford  a 


.;The  multipliers  1 

correction  which  takes  into  account,  as  a  first  aj.-proxi-  ; 
mat ion,  the  interpenetration  of  these  zones;  the  quantity: 

affords  a  correction  which  takes  into  ■ 
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aocountf  as  a  first  approximation,  the  intersection  of 
'  these  zones  which  arise  from  centers  of  various  types 
'(distributed'  with  a  uniform  density  and  in  forming). 

In  the  case  of  n  *  0  or  a  «  0,  we  get 


v(t) 


-  -f -  |(f  )!i 


(6) 

■:  In  the  first  case  the  initial  growth  of  the  product  is 
proportional  to  t^,  in  the  second  case  to  t*. 


B.  Ihe  formation  of  the  front. 

For  large  times,  when  t  ^t^^  and  t^  ,  the  evalua-  : 

.tion  of  the  integral  (2)  may  be  effected  by  using  the  so-  , 
■called  "pass  method",  which  takes  advantage  of -the  fact  ^ 
'that  the  function  under  the  integral  sign  is  appreciably  ■ 
different  from  zero  on3.y  in  a  small  neighborhood  of  the 
point  z  *  1.  Thus,  .we  insert  in  the  expression  under  the 

integral  sign,  1  -  ^  2(1  -  z),  1  -  +  2z^  5  3(l~z;)^ 

and  we  change  the  lower  integration.  limit  for  an  infinite-r 
;ly  remote  point.  After  some  calculations,  we  obtain: 


v(t)  ^  it  <1  - 


.^/gi 

^2  %/t 

dz7exp(■2^f) 

50 


(2)  ; 

The  integral  in  the  last  expression  is  equal,  up  to  the 
limit  of  a  . constant,  to  the  function  of  Kramp.  As  expres¬ 
sion  (2)  is  interesting  only  for  large  values  of  the  argu¬ 
ment  ,  so  using  an  asymptotic  expression  for  the  Kramp 
function,  we  obtain: 

1  -  +  1.5(— a.)  - 

From  the. latter  it  is  seen  that  as  time  goes  on,  the  ex-  ; 
pression  in  the  brackets  (curly)  tends  to  unity,  and 
v(t)“^it,  i«.,  a  flat  front  is  established,  which  spreads 
with  ..a  uniform  linear  yelocity  /t  .  In  the  case  when  „  „  : 


7 


and,  in  tbe  particular  case  v^hen  =  0,  or  a  «  0 


w(t)  =  X  — ~S»a»^  1  w(t;  =<  Air  w  ,  (i^) 

4  Aio 


As  is  seen,  for  la:c'?,e  tirr^es,  the  average  front  in¬ 
creases  v;ith  a  velocity  slis'htiy  lera-er  than  X.  This  is 
caused  by  the  fact  tbst  the  surl’acs  of  toe  real  front,  to 
which  the  transfonsation  propagates  perpendicularly,  is 
slightly  larger  than  the  area  of  the  final  plane  front. 

From  formulas  (10)  and  (^)  5  snd  equally  as  well 
from  figure  5,  it  is  seen  that  the  velocity  w(t)  should 
attain  a  aaxiitum  value  for  some  value  of  t. 

p.  Complete  form  of  the  curves  (2)  and  (5)« 

Let  us  consider  the  complete  form  of  the _ curves 
(2)  and  (3)  for  two  particular  cases:  for  n  =  0  and 
a  =  0.  In  this  investigation  it  is  convenient  to  change 
over  to  dimensionless  variables.  Thus,  v.'e  shall  obtain 
universal  curves,  Itoxu  v.'hich  all  experimental  curves  will 
be  obtained  for  difiersnt  a.,  and  a  by  the  mere  change 

of  scale.  In  the  capacit?;  of  a  time  interval,  we  shall 
take  the  intervals  introduced  above  tgc  and  t^^,  during 

which  time  the  front  of  the  propagation  of  the  transforma¬ 
tion  establishes  itself.  The  natural  scale  for  the  veloc¬ 
ity  w  is,  of  cour-se,  X*  Here  also,  the  natural  scale  for;, 
the  volume  will  be  equal  to  /tv,  or  Xt^*  the  voitaie 

transforoing  per  anit  ares  ia  the  time  t«<.  aad.  tj.,  for  the  frontal 
propagation  of  transformation  with  the  linear  velocity  . 

Thus,  for  n  =0,  we  introduce  the  dimensionless 


variables 


t-v4 


where  the  expressions  (2)  and  (^)  in  these  variables  are; 
f(^)  -  \e2cp£“t^(i  -  z)^(l  2z)_7ds/ 

1  (  0  ) 


1  —  ~  z)‘"(l  +  23)_7dz  + 


(15) 


hC. 

0 


+  2s)exp/-2^1“  z)^l  +  2z27dz. 
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Tl3.e  limit;  form  of  expression  (1^)  v/itti  dimension¬ 
less  variables  will  be;  “ 


"t  it.' 

f(r) 

1  M  1^)1  ‘V(T) 

«  ||2^)  (16) 

fir")  ^  t(i 

.  ■‘T/ 

“■^1^/12 /fi)i 

J!. 

V 

w(:2r)  i  + 

VlT/^8  t~  2  ,  (12) 

■  The  result  of  th.8  evaluation  of  vhe  integral  (15) 
for  all  values  is  sbcfe®  in  table  1;  the  corresponding 
curves  are  given  in  figures  4  and  5-  . 

The  velocitv  aiaoiimuia  lies  at  "T  ~  1-341  and  equals 

-■.!  O  T  X7.1  ■■  fii 


For  a  «  0»  introdxicing:  a  dimensionless  transforma¬ 


tion 


t 


rcnjl'^ 


n 


_v 


n 


V/ 


v-i  rrn^  ;  W  «  > 


(18) 


we  obtain 


g 


f  }  '■ 

1 

s#)  »  “  jexp/-^^( 

1  -  Z^)_7dZ  ( 

(  O' 

•  ( 

7 

w(^)  ■»  i  ~  ^  exp/-$^(l  - 


(19) 


-  z^)e-xp^^* 

2  ^(^) 
•(i  -  z  )7dz 


The  limit  values  of  the  formulas  (19)  end  (^) : 

P^ll  g(^)  »|^1  -  |^^)j  ^2C^(1  -  1^).  (21) 

>>1;  ^  ^  ^ 

Results  of  the  calculations  of  the  integrals  (2^)  and  (^) 
are  given  in  table  2.  In  this  case  the  maximum  velocity 
lies'^at  «  1.500  and  is  equal  to  W  *  1.285. 

Tn 


TABLE  1  TABLE  2 


E.  Position  and  magnitude  of  the  velocity  maximum 
for  the  general  case* 

We  are  investigating  the  change  of  the  velocity  as 
a  function  of  time  in 'two  particular  cases,  have  come  to 
the  conclusion  that  the  velocity  has  a  maximum,  and  have 

found  its  maximum.  , 

It  is  int^restirac j  now,  to  find  out  wbei^e  tnis  max— 
imum  velocity  shall  be  found  in  the  general  case,  and 
what  its  magnitude  will  be.  For  this,  we  will  present 
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the  integral  and  its  limit  values  (10)  and  (12)  in 
■the  folloijfing 

1  ‘  1  '  ■  '  '  ,  ■  '  ; 

i  ¥(r)  ='  i  -  -  z“)  ~g^(i  -  z)^(i  +  22)J?'dz  +; 


2^)  5^5(1—  2)^(1  +  22)^7exp^=1?|(i  -  z^) 


-  2)^(1  +  2z)_7dz;  (2^)  : 

at  r-^^i  . 

w(^)  «  2'^rii  “  |§1?)  +  “  M'f"'  (^i) 

at  i 


WCt;)  »  1  +3^ 


;  where  ii  a  dimensionless  pkraaetsr 

•  ^  .  y  '  x  -lip "" 

■  .  ■  6'  /  ?r2v^r  tir 

^  -  'f/  9  ~-V“  *  t' 

.  i  y  a"-  V 


(^)i 


;  defining  the  relative  density  of  center  of  the  propagation 
of  transfoiiaation  (distributed  with  a  uniform  density  with 
respect  to  formation  in  time).  It  is  easily  seen  that  tha 
parameter  ^  connects  the  dimensionless  parameters  intro-  : 
duced  above  ir  ftnd  ,  and  especially  , 

.  tDIie  result  of  the  calculation  X,  a^T  (?/)  and 

presented  in  table  3>. 

From  table  3  it  is  seen  that,  if  in  addition  to 
:  forming'  centers  of  transformation  with  constant  probabili- 
ity  of  incidence  a,  there  are  also  transformation  centers, 

;  distributed  with  a  uniform  density  n  ,  then  the  position  : 
i  of  the  maximum  velocity  tends  towards  smaller  values  of 
argument  2*  j  and  its  value  becomes  slightly  less.  If  the 
:  quantity  a  is  small  in  comparison  v/ith  n  ,  then  the  veloc- 
,  ity  maximum  is  determined  by  values  of  argument  'i'- close 
;to  zero,  and  its  magnitude  is  slightly  different  from 
the  value  vr  *  1.285* 

IB 

In  the  converse  situation,  if  to  the  transformation 
■centers  distributed  with  constant  density  n  are  added 
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X i -  "t.  1  r  j  XT*  X  ?i?  ’O s’‘*3i.S7)*'  ^O*  ^  XAS^‘X.ri^-. 

table'^,  ond  t-ho  deperjder'ce  ]i^=  TW,  czxlj  in  this  esse  the 

■oarsmeter  sho’jld  I/>/  .  v.  -v 

'iTi  fS-injrs  4  the 'bold  line  rer'r-esei;t5  C’jrve  C xt j  ^ 
btA  the  dashed  line  is  the  oui'-ve  rbteined  froa  Integral 
■<2) ,  eicpx-essad  in  di&-:sr-6ior,.leos  variables: 

f  (T)  ^t/i  “  \  -*  z*^)  -r-'Tl  “  z)^* 

/  .  i  .  ~  )  _ _ _ 


/.  v-i-  •'  n.!.  jL 


(27) 


for  c'  4-1 

i'C  r  ) 


-  ^-.•i  7 


/  %  yi  t.  c  >7 


2,22  15  ^3 


5^t^I.7(2§) 


for  7'  ^“1 


frr)  »  r/i 


-  (1 


2  ry/^ 


27  _ _  iTCPQl 


and  cslcnlatsd  for  'n  ^  1  *  In  t'lgn.re  y  tne  bole  Ixne  is 
the  cvir/e  of  the  velocity  frniri  (V?),  obtained  r-Oirv 
for  the  values  of  the  paraaeler  O.b  (first  assr-ed 
curve)  and  7“  ~  >  (the  sacoxji  dssbea  curve). 

5.  The  €3{:ria.TiS;C::tai  detennlnatioa  of  the 
liiagrTiTudes  oT^T“i'\^  ■)  ana  a. 

The  entities  and  n  can  be  easily  deternunsd 
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experimental l.y.  i:or  an  observed  velocity  of  propagation 
of ■  transfoririaticn ,  after  an  -appreciable  amount  of  time  has 
elapsed  since  the  baginnir.g  o.f  tbs  process,  when  already 
to  all  intents  Cincl  pui-pyosea  a  plane*  front  has  been  estab¬ 
lished,  we  will  have 
'  -  Z  -  li®-  w(t)  . 

'  ,  t~->CI>  - 

The  velocity  in  the  initial  period  of  the  process  enables 
us  to  determine  ,  knowing  vi , 


1  w{t)  b 


n  a  - ^  iim  . 

^  2T\P  t '-■^0 

In -the  case  v/hen  -  0,  similarly  we  obtain  a; 

o 
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